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Abstract 

A peculiar property of complex scattering potentials is the appearance of spectral singu- 
larities. These are energy eigenvalues for certain scattering states that similarly to resonance 
states have infinite reflection and transmission coefficients. This property reveals an interest- 
ing resonance effect with possible applications in waveguide physics. We study the spectral 
singularities of a complex barrier potential and explore their application in designing a waveg- 
uide that functions as a resonator. We show that for the easily accessible sizes of the waveguide 
and its gain region, we can realize the spectral singularity-related resonance phenomenon at 
almost every wavelength within the visible spectrum or outside it. 
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1 Introduction 

Recently it has been noticed that by defining the Hilbert space of a quantum system using a 
nonstandard inner product one can generate a unitary time-evolution by a Hamiltonian operator 
that is manifestly non-Hermitian in the standard L^-inner product [1]. A typical example is if = 
— ^ +ix^. It turns out that one can actually describe the corresponding physical system using the 
standard Hilbert space and a Hermitian Hamiltonian which is much more complicated and highly 
nonlocal [2] . The non-Hermitian Hamiltonians with the above property must have a real spectrum 
and a complete set of eigenfunctions so that every wave function can be expanded in terms of these. 

*E-niail address: aniostafazadeh@ku.cdu.tr (corresponding author) 
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There are two main mechanisms that are responsible for the incompleteness of the eigenf unctions of 
a non-Hermitian Hamiltonian. These are respectively associated with the presence of exceptional 
points [3] and spectral singularities 

An exceptional point is a point in the space of parameters M of the Hamiltonian operator H 
where an eigenvalue of H becomes defective, i.e., if the parameters change continuously along a 
path in M that passes through an exceptional point, two or more of the eigenvalues of H together 
with their eigenvectors coalesce. Exceptional points have found various physical realizations and 
applications [5]. They have also been the subject of experimental studies [6]. In contrast, spectral 
singularities are certain points of the continuous spectrum of non-Hermitian scattering Hamiltoni- 
ans whose presence ruins the completeness of the eigenf unctions, although for each point of the 
spectrum (including the spectral singularity) there corresponds two linearly independent eigenfunc- 
tions The physical meaning and importance of this strange mathematical phenomenon have 
been obscure until very recently [8]. 

The appearance of spectral singularities as an obstruction for defining a unitary quantum system 
using a complex potential was initially noted by Samsonov [9] who, following the pioneering work of 
Naimark |11[10], only considered systems defined on the half-line. A simple example of an exactly 
solvable model defined on the whole real line that can support a spectral singularity is the one 
given by a delta-function potential, v{x) = ^S{x), with a complex coupling constant 3. As shown 
in [H] this model has a spectral singularity provided that 3 is purely imaginary. The emergence 
and the mechanism by which a spectral singularity spoils the completeness of the eigenfunctions of 
a scattering non-Hermitian Hamiltonian are studied in [7]. 

In Ref. [Sj we have provided a physical interpretation for the spectral singularities that identifies 
them with the energies of scattering states having infinite reflection and transmission coefficients. 
Because this is a characteristic property of resonance states, spectral singularities correspond to 
the resonance states having a real energy and hence a zero width. The existence of such states 
reveals a new type of resonance effect with potential applications in various areas of physics. To 
demonstrate how one can realize this effect, we have explored, in Ref. [8], the occurrence of spectral 
singularities for the PT-symmetric imaginary barrier potential: 



where ( and a are real parameters, a is positive, and sgn(x) stands for the sign of x n The 
potential ([T]) can be used to model the propagation of certain transverse electric fields in a planar 
slab waveguide with an adjacent pair of loss and gain regions [I2] as shown in Figure [H In [8], 
we have determined the values of physical parameters of this system that correspond to a spectral 
singularity. It turns out that for fixed values of the parameters entering the permittivity, spectral 

^Here we only consider Hamiltonians of the form H = — 'Si^'^'^i^) where x is real and v{x) is a complex scattering 
potential. 

^ sgn(a;) :— x/\x\ if a; 7^ 0, and sgn(O) := 0. 
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Figure 1: (Color online) Cross section of a rectangular waveguide with gain (+) and loss (-) regions 
in the x-z plane. Arrows labeled by /, R, and T represent the incident, reflected, and transmitted 
waves. 

singularities occur for a discrete set of values for the length of the loss and gain regions 2a and 
the height of the waveguide 2/3. For the fundamental mode of the standing wave formed in the 
transverse direction (along the x-axis), the largest allowed values of a and /3 are respectively 
1004.17 nm and 62.0464 nm. Explicit calculation shows that for these value of a and (3, the 
reflection and transmission coefficients (at the resonance frequency) exceed 10^^. Note however 
that this requires maintaining this value of a for the length of the loss and gain regions to an 
extremely high degree of precision (10~^ Angstroms). A 0.1% deviation from the above values 
of a and /3 leads to a decrease in the reflection and transmission coefficients by a factor of 10''. 
Furthermore, the experimental realization of spectral singularities that is proposed in [8] involves 
restricting the frequency of propagating wave to the resonance frequency of the active region and 
demands that the imaginary part of the permittivity for the loss and gain differ only by a sign. 
In the present article, we study the spectral singularities of the complex barrier potential 



where 3 is an arbitrary complex coupling constant. Our main motivation is that the proposed 
experimental setup that would confirm the resonance phenomenon associated with the spectral 
singularities of (|2]) is free of most of the above-mentioned restrictions. This is essentially because of 
the additional degree of freedom, namely the real part of 3 in ([2]). Another reason for considering the 
potential ([2]) is that, similarly to ([T]), the mathematical problem of locating the spectral singularities 
of ([2]) admits an exact and essentially analytic solution. This is very convenient, for it enables us 
to avoid dealing with the subtleties of the numerical treatments of the problem. 

A most remarkable result of our investigation is that for the experimentally easily accessible 
dimensions of the waveguide associated with ([2]), we can produce the spectral singularity- related 
resonance effect at almost every wavelength within the visible spectrum or outside it. This is of 
great practical significance, for it may form the basis of a new mechanism for producing laser beams 
of desired wave length. 

This article is organized as follows. In Section 2, we provide a brief review of the results of 
[3 [8] on the calculation of spectral singularities and their relation to the refiection and transmission 




for |x| < a, 
for |x| > a, 



(2) 
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coefficients. In Section 3, we solve the problem of locating the spectral singularities of the complex 
barrier potential. In Section 4, we describe a waveguide modeled using this potential and explore 
the possibility of the detection of the resonance effect associated with the spectral singularities. In 
Section 5, we present our concluding remarks. 



2 Spectral Singularities 

Consider the Hamiltonian operator H = —-^ + v{x) for a complex scattering potential v{x) with 



x G M. Suppose that v{x) tends to zero as x — )■ ±oo so rapidly that the integral j _^{l + \x\)\v{x)\dx 
converges. Then the spectrum of H has a real continuous part. Suppose for simplicity that the 
spectrum coincides with [0, oo). Then the eigenvalue equation Hip = Eijj yields the following 
asymptotic expressions for the (generalized) eigenfunctions of H. 



iplix) A^e'^'' + Ble-'^"", for x ±oo. (3) 



Here g is a degeneracy label taking values 1 and 2, k := ^/E, and A^., B\. are possibly fc-dependent 
complex coefficients that are related by the so-called transfer matrix M = (Mjj) according to 



Bi 




(4) 



The transfer matrix encodes all the necessary information about the scattering properties of the 
system. It is easy to show that it has a unit determinant [7]. 

Next, consider the Jost solutions ipk± of the eigenvalue equation Hip = k'^ip that by definition 
satisfy 

ipk±{x) -f e^''^'' for X ±oo. (5) 

Denoting the coefficients Aj_ and B^ for the Jost solutions ipk± by and B^ and using of ([3]) - 
dS]), we find AX = BZ = 1, Az = B+ = 0, At = B^ = M22, A; = M12, and Bl = -M21, [7]. 
These relations identify ipk± with the left- and right-going scattering solutions [13]. As a result, 
the left and right transmission T'''' and reflections i?''^ amplitudes are given by [8] 

M22 M22 M22 ^ ' 

Spectral singularities of H are eigenvalues E = k"^ for which ipk+ and ipk- become linearly- 
dependent [H [m |T5]. They are given by the real zeros of M22 where both the transmission and 
reflections coefficients diverge [8]. Because the latter is a characteristic feature of resonances [16] . 
we associate spectral singularities with a peculiar type of resonances which satisfy the eigenvalue 
equation Hip = k'^ip for a real k"^. These are resonances having a zero width, [8]. 
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3 Complex Barrier Potential 



Consider the complex barrier potential (j2]). The solution of the time- independent Schrodinger 
equation, Hip = k^ip, for this potential is elementary. We can use this solution, to determine the 
entries of the transfer matrix M: 

Mil = -r^ (7) 



Mi2 = -M21 



Aw 

i{w^ — 1) sin{2akw) 
2^ 



where w := a/1 — 3//^^ §, and for all x ^ Il^i 

f{w, x) ■■= e-''^"(l + wf - e''^"'(l - wf. (10) 

According to ([9]), spectral singularities are associated with real values of k for which 

f{w, ak) = /(v/l-3A^ = 0. (11) 

This is a complex transcendental equation involving four real variables a, k, Re(3) and Im(3), where 
"Re" and "Im" stand for the real and imaginary part of their argument. A direct numerical solution 
of ( ITT]) would amount to identifying two of the variables as parameters (say Re(3) and Im(3)) and 
solving ( fTTl) for the remaining two variables (i.e., a and k). This requires selecting a range of values 
for the parameters and performing a numerical solution of (fTTj) for all the values of the parameters 
within this range. This is not a straightforward procedure, for in principle there is a real solution 
only for a one-dimensional subset of the two-dimensional parameter space. Numerical treatments 
may miss part of this one-dimensional subset and could lead to various types of errors @. In the 
following, we shall avoid these difficulties by offering an exact and essentially analytic solution of 
(|TT]) . Specifically, we shall reduce (|TT]) to a single real transcendental equation involving two real 
variables and a discrete label. As we will see this provides a much more efficient and precise method 
of locating spectral singularities of the potential (E]). 

We begin our analysis by recalling that in light of ( ITOj) . we can express (fTT]) in the form 

e^"^"'(l - w) = ±€-^"'="'(1 + w). (12) 

This is equivalent to 

cos(2aA;v/l -3^2) = ± i — - 1 J . (13) 



■^We use -^/l — 3/fc2 to denote the principal value of (1 — j/fc^)^/^ whose argument belongs to [0,7r). 
*For example typical numerical shemes will produce k values with a small imaginary part. These do not correspond 
to scattering states and must therefore be rejected. 
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We can reduce (fT3|) into a pair of real equations. To this end, we first introduce 

Re(3) ^ Im(3) 



cr 



P 

Q 
r 



a :- 



A;2 ' 



y 



P 



at 



/cV2|l-p|(vV+T-l) sgn(y), 



: ak^J 2\l- p\{.Jy^ + 1 + 1). 
In light of (fT5l) . ( IT61) . and /c > 0, we have 

sgn(g) = sgn(?/) and r > 0. 



(14) 

(15) 
(16) 

(17) 



Next, by inserting i/k"^ = p + ia into ( IT3l) and using the elementary properties of the trigonometric 
and hyperbolic functions, we express the real and imaginary parts of this equation as 



cosr coshg 
sinr sinhg 



(i-p)V + P^ ^ 

^ 2(1 -p)y 
+ (l-p)V + p'' 



(18) 
(19) 



Here either the top or the bottom sign should be taken in both of the equations. 

We can easily eliminate r in f fTSj) and f|T9l) . We do this by solving for cosr and sinr and 
imposing the identity sin^ r + cos^ r = 1. Expressing cosh^ q that appears in the resulting equation 
as 1 + sinh^ q, we find after some lucky cancelations 

4(1 -p)V 



sinh q 



4(1 - p) 



sinh q 



,(i-p)V + p2 

Solving this equation for q and using f|T7|) . we obtain 

q = sgn{y)Q{p,y), 

where 

Qip, y) ■= sinh"^ 



[(1 _p)2y2 + ^212 



0. 



(20) 



2(v/y2 + i|i-p| + i-p) 



(l-p)V + p' 

Next, we substitute ( 120|) in ( ITSl) and solve for r. Employing the identity cosh(sinh~^(x)) 



(21) 



VT+~a? to simplify the resulting expression, we then find 

r = Rnip,y), 

where n is an integer, e is a sign, and 



RniP^y) ■= ™ ± cos ^ 
Three remarks are in order. 



L-ll-ply/^^ ' 

v/(i-p)V + P^ 



(22) 



(23) 
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1. Because for real values of p and y, we have 



|l-p|v/y=^ + l + p-l>0, (24) 
(1 - p)Y + = (l - 1 1 - p| ^2/2 + 1 j ' + 2 (|1 - p| ^2/2 + 1 + p - l) , (25) 

the argument of cos~^ in f l23|) lies in the interval [—1, 1]. This implies that the second term 
on the right-hand side of (l23i) takes values in [0, vr], and R^ip, y) is always real. Furthermore, 
sgn(i?^(p, y)) = sgn(n) for n 7^ 0, and sgn(i?^(p, y)) = ± except for the case p = y = 
which corresponds to a real barrier potential. Therefore, to ensure that r > 0, we must 
demand that either n > and e be arbitrary or n = and e = +. In the following we shall 
confine our attention to these two cases. 

2. Every r and q that satisfy f lTSj) and ( 1T9|) must also satisfy (120|) and f l22|) . but the converse is 
not true. It turns out that depending on the choice of the top or bottom sign in ( ITSl) and 
(fT9|) and range of values of p and y one must further restrict the values of n and e in ( l22l) . 



3. Equations (1181) and ( 1191) are not equivalent to ( 1131) unless we impose ( 1151) and ( 1161) . First, we 
impose ( 1T6|) . This allows us to express ak in terms of r, p and In view of (|22|) . we then find 

afc = G^(p, y) := (26) 



2|l-p|(V^^+l) 

Unlike ( |T6i) that determines ak in terms of p and y, (|T5l) gives rise to an additional constraint 
on the possible values of p and y. We will next derive an explicit form of this constraint. 

We can use (|T5|) and (|T6|) to obtain 

^ = ^\ r-o T — 7 ^Sn(y) . (27) 

If we insert ( !22l) in this equation and use (fT7j) , we arrive at 

q= sgn(y)g^(p,|/) (28) 

where 

Qnip.y) ■.= K{p.y)\h^^^^. (29) 

Recall that here n is a nonnegative integer, and for n = we have e = +. 

Enforcing ( l20l) and ( 128|) yields Q{p,y) = Qn{p,y)- This gives an infinite sequence of equations 
that determine the locus of points in the p-y plane for which the complex barrier potential has 
a spectral singularity. These form an infinite sequence of curves. Clearly, Q{p,y) = QniPyV) is 
equivalent to F^{p,y) = where 

Kip, y) ■■= \ [sinh^ Q(p, y) - sinh^ Q^(p, y)] . (30) 
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In view of fl2T|) . f l29|) . and fl30|) we have the following explicit expression for F^. 
|l-p|v/^^+l-p 




7m ± cos ^ 



v/(i-p)V + p' 



(31) 



In the p-a plane, which coincides with the complex 3/A;^-plane, the spectral singularities are located 
on the curves 

(32) 



C::={(p,a)GM2 |f,^(p,^)=0}. 



As we will see below, not all the points on these curves turn out to correspond to a spectral 
singularity. 

Examining the graphs of C^, see Figure [21 we find that these curves all lie to the left of the line: 
p = 1 in the p-a plane. The curve Cq intersects this line at (p = 1, a = 0) but this point does not 
correspond to a spectral singularity, because for o" = the potential is real. These observations 
allow us to confine our attention to the case p < 1. A consequence of this relation and fll7p is 



sgn(g) = sgn(t/) = sgn(o-) = sgn(Im(3)). 



(33) 



In order to determine the choices of n and e for which includes points corresponding to 
spectral singularities, we proceed as follows. First, we recall that the spectral singularities are give 
by ffTT]) . Next, we use the identities w = y/1 — p — ia = ^1 — p — i(l — p)y, and Eqs. ( fTOl) and 
to express the left-hand side of ffTTl) in terms of p and y. This yields a sequence of equations 
of the form F^{p, y) = 0, where 



and Gt is introduced in 



-2iG^ {p,y)^l-p-i{l-p)y 



l + p-iil-p)y 



_ 2iGl,{p,y)^l-p~i{l~p)y 



1 - v^l - p - i(l - p)y 
The spectral singularities lie on the curves 



(p,a)GR^ |F^(p,^)| = 



(34) 



It turns out that |-f^(p, -^)\ = only for e = — and n > 0. Furthermore, as shown in Figure^ the 
graph of C~ coincides with a part of the graph of that lies above the p-axis in the p-a plane. 
That is 

(35) 



C- = ^{p,a)eM.' F„-(p,^) = 0, a>0 

This provides a highly nontrivial check on the validity of our calculations. It also reveals the 
interesting fact that spectral singularities can appear only for the cases that the coupling constant 
3 has a positive imaginary part. In Section 4, we will provide a simple physical justification for this 
fact. 
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Figure 2: (Color online) Graphs of the curves C~ (the dashed black curves) and C~ (the solid 
red curves) in the p-a plane for n = 1,2,3. C~ are symmetric about the p-axis. coincides 
with for cr > 0. This shows the lack of spectral singularities for a < 0. The dashed black line 
corresponds to p = 1. The solid gray lines are the coordinate axes. 
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Figure 3: (Color online) Graphs of the curves C~ in the p-a for 1 < n < 20. The unmarked curves 
correspond to n = 5, 6, 7, ■ ■ ■ ,20 from top to bottom, respectively. Cf has an asymptote given by 
p ^ 0.667. The dashed line corresponds to p = 1. The solid gray lines are the coordinate axes. 

Another outcome of examining the curves C~ is that, in comparison to f l35|) . the use of flM|) 
in ploting C~ leads to substantially larger numerical errors. Therefore, in practice, it is more 
appropriate to employ f l35|) to locate the spectral singularities of the potential ([2]) graphically. For 
this reason, in preparing Figure [3l we have used (l35l) to plot C~ for 1 < n < 20. This figure reveals 
the following facts about the curves of spectral singularities. 

• Ci has a vertical asymptote given by p 0.667. 

• For all n > 1, C~ is a decreasing curve. 

• For n > 2 the graph of C~ lies above that of Cn^i- 

• As ?T, — )■ oo, C~ tends to the p-axis, §. 

4 Resonance Effect Associated with Spectral Singularities 

In this section, we describe an electromagnetic waveguide that can be used to realize the resonance 
effect associated with the spectral singularities of the complex barrier potential 02]). 

^This can be easily shown using (HU)). 
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Figure 4: (Color online) A rectangular waveguide aligned along the z-axis. The dashed region is 
filled with an atomic gas. i, j, k are the unit vectors along the x-, y-, and 2-axes, respectively. 

Let a,/3,7 be positive parameters having the dimension of length. Consider a rectangular 
waveguide with perfectly conducting walls that guides transverse electric (TE) waves along the 
2;-axis in the region defined by \x\ < (3 and \y\ < 7. Suppose that the region \z\ < a inside the 
waveguide is filled with an atomic gas so that the relative permittivity inside the waveguide is given 

by mm 



2i5uj 



for \z\ < a, 
for Ul > a, 



(36) 



where u,Up,uo and S are respectively the frequency of the wave, the plasma frequency, the resonance 
frequency, and the damping constant. Figure H] provides a schematic illustration of the waveguide. 

Let ^ := u/c be the wave number, '■= i^m/{2P) for every positive integer m, and (j){z) be a 
solution of 

0"(;.) + [^h{z) - ^l]<P{z) = 0, (37) 

with a continuous derivative. Then it is an easy exercise to show that the following TE wave is a 
solution of the Maxwell equations with appropriate boundary conditions for the above waveguide 



E{r,t) 
B{f,t) 



iu sin[j^m(a; + /?)] (/)(z)^j| , 

i[J^(x + /5)] (l)'{z)i - iim. cos[^{x + /?)] (I){z)k^ } 



(38) 
(39) 



Here £^ is a real constant and i, j, k are respectively the unit vectors along the x-, y-, and 2-axes§. 
Eq. (137|) coincides with the eigenvalue equation. 



dz^ 



+ v{z) 



= fcV, 



^Note that we use physicists' convention of expressing the time-harmonic part of the wave as e"*"^* as opposed 
to engineers' convention that uses e'"*. This leads to a change of sign in the imaginary part of permittivity (and 
complex refractive index.) 
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for the complex barrier potential ([2]) , if we replace x with ^ in ([2]) and make the following identifi- 
cations. 



..yi^.^^i-^, ,40) 



■ i 



These imply 

P 



a 



(43) 



For a TE wave given by (j38l) and (j39ll that propagates in the waveguide, we have ^ > so that 
k is real and positive. Therefore, we can use the results of Section 3 to determine the reflection and 
transmission amplitudes of the wave. In particular, if we can arrange the physical parameters of the 
system namely u},Uo,u}p,6,a and /3 so that k"^ is a spectral singularity, reflection and transmission 
coefficients diverge and the waveguide functions as a resonator. 

For the system we described above Up and S are positive. In view of (l42l) this implies that 
Im(3) < 0. Therefore, as we showed in Section 3 (Figure [2]), there is no spectral singularities. 
This is actually to be expected, because for the physical system that we consider the presence of a 
spectral singularity means the emission of highly amplified electromagnetic waves. This is clearly 
forbidden by the law of conservation of energy. It is quite remarkable that the absence of spectral 
singularities for the barrier potentials with Im(3) < 0, which is a purely mathematical result, can 
be interpreted as a direct consequence of conservation of energy in classical electrodynamics. 

We can avoid the above conflict with energy conservation, if we consider the situation that the 
sign of Up is reversed. In this case the atomic gas confined in the region < a inside the waveguide 
acts as a gain medium [18] . We can achieve this by shining a laser beam along the ?/-axis to excite 
the resonant atoms and induce a population inversion 0. This was also the mechanism we employed 
in [8] to realize the resonance effect related with the spectral singularities of the PT-symmetric 
barrier potential ([1]). The main advantage of the system we consider here is that, for the cases that 
uq, Up, and 6 are fixed, the spectral singularities occur along an infinite set of curves in the relevant 
parameter space, whereas for the system considered in [8] they occur at an infinite set of isolated 
points. 

Another practical shortcoming of the system studied in [8] is that it relied on the assumption 
that the frequency of the propagating wave u is equal to the resonance frequency of the gain medium 



is proportional to the difference between the number density of atoms in the ground and excited states which 
becomes negative due to population inversion [18) . 
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ujq. Here we do not need to make such an assumption. This provides an additional flexibihty in 
our attempt to adjust the parameters of the system so that the frequency of the propagating wave 
coincides with that of a spectral singularity. To demonstrate how we can do this, we first choose 
some typical values for oJoyUJp and 6 that determine the characteristics of the gain medium. For 
example, suppose that 

huo = 5 eV, h'^ul = -0.04 eV^, M = 1.25 eV, (44) 

Then UQ,ujp and 6, p and a become functions of Q and u. For each choice of Q, we can select a 
practically attainable range of frequencies u and plot the parametric curve defined by (143|) in 
the p-a plane for this range of values of u. The points of the intersection of and C~ correspond 
to the values of u for which k'^ is a spectral singularity provided that (|26l) is also satisfied. The 
latter condition fixes the value of a. The u and a values obtained in this way depend on Q, the 
label n, and another discrete label that counts the number of times intersects C*^ as one 
decreases the value of p starting from its largest allowed value, namely 1. 

To determine the precise values of u and a for a given choice of Q and n, we plot the curves 
and C~ in the p-a plane as in Figure O and select an intersection point, i.e., fix in- We then 
determine the coordinates (p^,crv,) of this point graphically §. Next, we numerically solve for u in 

p = pi,, cr = (T^, (45) 

where the left-hand side of these equations are given by ( l43l) . A consistency check on this calculation 
is that both of these equations must give the same value for u. Having obtained we compute k 
using (l40l) . substitute this value of k and (145|) in (!26|) . and solve for a in the resulting equation. 

Figure E] shows the plots of the curves and C~ for various values of Vt and n. For every fi, the 
curve intersects each of the curves C~ at least once and at most three times, i.e., 1 < < 3. 
As one adopts sufficiently large values for /3/m (sufficiently small values for (]), the number of 
intersection points increases. This number equals to the number of different frequencies u at which 
a spectral singularity occurs. Note, however, that each of these spectral singularities correspond to 
a different value of the parameter a. Furthermore, because for each value of f2, intersects C~ 
for all n, for every value of /3/m there is an infinity of choices for the u (and a) values that yield a 
spectral singularity. 

It turns out that for each value of /3/m, one can select n so that the related spectral singularities 
correspond to experimentally preferable values for a and uj. For example, for the principal mode 
(m = 1) of a waveguide of height 2/3 = 1 cm, choosing n = 10000 and in = 2, we find a spectral 
singularity of frequency u ^ 2.1555 ev/h (wave length: A ~ 575.20 nm) provided that we maintain 
an active region of length 2a ~ 2.8786 mm. For future reference we will use to label this 
spectral singularity. In fact, as we would expect, changing the value of n starting from such a large 

^Using Mathematica and without making any special efforts we could zoom in the neighborhood of the relevant 
intersection point and obtain and ct* up to six significant digits for various values of Jl, n, and in- 
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Figure 5: (Color online) Graphs of (thin black curves) and C~ (thick colored curves) for 
Uo,Up,6 as given by ( l44l) and Q = 1.0,2.0,2.5,3.0 eV/h that respectively correspond to f3/m = 
310.0, 155.0, 124.0, 103.3 nm. The top figure, shows the intersection points of and the C*" with 
n = 1,2,3,4,5. These correspond to in = 1- The bottom figure shows the intersection points of 
and C- with n = 600, 800, 1000, 2000, 3000. For these points = 1, 2 or 3. 
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Figure 6: (Color online) Location of the spectral singularities in the 2a-\ plane for 2/3 /m = 1 cm 
and 10000 < n < 11000. The numerical values of 2a (length of the gain region) and A (wavelength 
of the emitted wave) are given in mm and nm, respectively. The displayed dots represent the 
spectral singularities with n = 10000, 10100, 10200, ■ ■ ■ , 11000 from the left to the right. The curve 
segments joining adjacent dots represent 99 isolated points each corresponding to a different spectral 
singularity. The top and bottom figures show the spectral singularities with = 2 and in = 3, 
respectively. The spectral singularity is also displayed. 



initial value, we find a large number of closely spaced spectral singularities all belonging to the 
experimentally attainable range of values of a and u. 

Figure [6] shows the location of the spectral singularities for 2/3/m = 1 cm that we obtained for 
10000 < n < 11000 and in = 2, 3. For these choices of n, the spectral singularities associated with 

= 1 correspond to extremely large (experimentally unattainable) values for 2a. As suggested by 
Figure |6l using a waveguide of height 2/3 = 1 cm and adjusting the length 2a of the gain region 
(in the millimeter range), we can generate spectral singularities at almost every point of the visible 
spectrum or outside it. This is also supported by a detailed examination of other choices for the 
values of 2/3 and n that we do not fully report here, except for a few illustrative cases that we 
include in Tables [1] and [2l 

As shown in Figure |6] and Tables [1] and [21 for in = 2 the wavelength of spectral singularity is 
an increasing function of the length 2a, while for = 3 it is a decreasing function of 2a. This 
shows the availability of a wide range of values of parameters of the system that lead to a spectral 
singularity. 

According to Table [H for n = 10000 and = 2, both of the choices 2(3 /m = 1 mm and 
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2/3/m = 1.0000 





A 


2a 




1 


1.6798 iim 


15.517 mm 


0.99919 - 6.1408 x lO'^i 


2 


614.03 nm 


3.2291 mm 


0.99910 - 2.1814 x lO'^^i 


3 


162.61 nm 


0.81531 mm 


1.00045 - 2.6078 x 1{)-H 



213 /m = 1.0000 mm 





A 


2a 




1 


1.9974 mm 


307.845 m 


0.99920 - 4.9699 x 10"*^ 


2 


575.20 nm 


2.8786 mm 


0.99908 - 2.4333 x 10~^i 


3 


162.85 nm 


0.81379 mm 


1.00045 - 2.6259 x lO^^i 



2/3/m = 1.0000 cm 



f 


A 


2a 




1 


1.9982 cm 


6.7894 km 


0.99920 - 4.968 x 10~^i 


2 


575.20 nm 


2.8786 mm 


0.99908 - 2.4333 x lO^^i 


3 


162.84 nm 


0.81379 mm 


1.00045 - 2.6259 x lO'^^i 



Table 1: The wave length A of the TE wave, and the length 2a and complex refractive index ^/e of 
the gain region corresponding to a spectral singularity with ijJp,uo, S given by (l44l) and n = 10000. 
The boldface figures are those of the spectral singularity i^. 











n 


A (nm) 


2a (mm) 






2000 


306.59 


0.30685 


0.99902 - 1.1437 x 10"^^ 


3000 


347.47 


0.52173 


0.99893 - 7.763 x 10"^^ 


4000 


382.28 


0.76534 


0.99895 - 5.8934 x lO'^^i 


5000 


415.09 


1.03877 


0.99897 - 4.757 x W^i 




in — 3 






n 


A (nm) 


2a (mm) 






2000 


220.78 


0.22059 


1.00055 - 1.1701 X lO-^i 


3000 


203.54 


0.30504 


1.00064 - 8.043 x lO'^i 


4000 


193.10 


0.38589 


1.00062 - 6.1592 x lO'^^i 


5000 


185.45 


0.46327 


1.00059 - 5.006 X lO'^i 



Table 2: The wave length A of the TE wave, and the length 2a and complex refractive index y/e 
of the gain region corresponding to a spectral singularity with ujp,ujQ,5 given by pil) . 2(3 /m = 
1.0000 cm, = 2, 3, and n = 2000, 3000, 4000, 5000. 
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2/3/m = 1 cm give rise to spectral singularities with the same values for the wave length A, the 
length parameter 2a, and the complex refractive index ^Je. Therefore, to observe this spectral 
singularity one does not need to fine tune the height of the waveguide. This is also true for 
n > 2000 and C = 3. 

Figure [7] shows the graphs of logiodT'P + l-Rp) as a function of u/us where Us is the frequency 
of the spectral singularity with n = 10000 §. This is represented by the peak at w/ws = 1 that 
corresponds to an amplification of the emitted electromagnetic energy by a factor of > 10^^. We 
can increase this factor, if we use more precise values for the u and 2a parameters associated with 
The other peaks represent the spectral singularities with n = 10001, 9999, 10002, 9998, ■ ■ ■ and 
= 2. These are lower than the central peak, because the corresponding spectral singularities 
have a slightly different values for u and 2a than those for that we have adopted in preparing 
this figure. The location of the peaks can be easily explained. Because the frequency of the spectral 
singularities with n = 10000 and = 2 decreases as a linear function of n, changing n by one unit, 
i.e., by a factor of 10~^, leads to a change of u by the same factor. This explains the occurrence 
of the peaks at the u/ug values that differ from unity by an integer multiple of 10"^. Another 
interesting observation is that the amplification effect persists even if we change u by as much as 
10% (See the last graph depicted in Figure [3) This is another indication of the suitability of the 
proposed setup for an experimental realization of the spectral singularity-related resonance effect. 

5 Conclusion 

Spectral singularities are the energies of a rather strange type of scattering states that similarly to 
resonance states have infinite reflection and transmission coefficients. Physically this means that 
if one arranges to scatter a plane wave of energy close to that of a spectral singularity, the system 
amplifies the wave enormously. Therefore, spectral singularities are related to a novel resonance 
phenomenon that awaits experimental verification. 

For a real potential, the unitarity of dynamics implies that the refiection and transmission 
coefficients add up to unity. Therefore, spectral singularities cannot arise for a real potential. In this 
article, we have explored the spectral singularities of a complex barrier potential and showed how 
this potential enters the description of the electromagnetic waves propagation in certain waveguides. 
If we adjust the physical parameters of the waveguide system so that the frequency of the incoming 
wave is close to that of a spectral singularity, the waveguide emits a substantially amplified wave 
of the same frequency from both ends. We expect this phenomenon to find useful applications in 
optics, particularly because we can adjust the resonance frequency to attain any desired value. 

A remarkable result of our investigation is that spectral singularities of the complex barrier 
potential occur only for certain values of the coupling constant with positive imaginary part. While 

^Because of the symmetry of the problem, and as we can easily verify by inserting ([5]) in ([5]) , the left and right 
reflection (and transmission) coefficients coincide; = =. R and T'' = T' =: T. 
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Figure 7: (Color online) Graphs of log]^Q(|Tp + |-Rp) as a function oi lo/uJs for three different ranges 
of values of oj/uJs- |-Rp and |Tp are the reflection and transmission coefficients (See endnote [27].) 
Us = 2.15548 ev/h is the frequency of the spectral singularity ir that occurs for Up,UQ,6 given 
by dSD, n = 10000, 4 = 2, 2/3 /m = 1.00000 cm, 2a = 2.878644 mm, A = 575.2046 nm, and 
y/e = 0.999081358 — 0.0002433302. The highest peak represents The other peaks correspond to 
the spectral singularities with n = 10001, 9999, 10002, 9998, ■ ■ • . 
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the mathematical reason for this phenomenon is obscure, it has a simple physical justification. The 
condition that the coupling constant has a positive imaginary part corresponds to the case that 
the waveguide includes a gain region. This is necessary for the existence of the resonance effect, 
because of the law of conservation of energy. 

One of the main goals of the present investigation was to examine the range of critical values of 
the physical parameters of the waveguide system that corresponded to the occurrence of spectral 
singularities. We showed that indeed these parameters spanned a very large spectrum of values. 
This makes the waveguide system we considered in this article a plausible candidate for observing 
the spectral singularity-related resonance effect. 
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